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After a short introduction on some general implications of recent cosmological observations, we discuss the 
properties of relic neutralinos in connection with present experimental strategies for detecting WIMPs. 



1. Introduction 

A large number of independent astronomical 
observations (on galactic halos, clusters of galax- 
ies, large scale structures) indicate that in our 
Universe the average density of matter, normal- 
ized to the critical density p c = 1.88 • 10~ 29 h 2 g • 
cm~ 3 (h is the Hubble constant in units of 100 
km • s" 1 -Mpc -1 ), is in the range 0.2 < Q m < 0.4 
(or equivalently, 0.05 < Vl m h 2 < 0.3) 0. The 
value Q m ~ 0.3 may also be derived by combin- 
ing CMB measurements 0| with data on high- 
redshift SNIa 0|, with the further result that a 
large value, SIa ~ 0.7, should be assigned to the 
cosmological-constant (or quintessence) contribu- 
tion. 

By comparing these results on f2 m with the 
contribution to matter provided by visible mat- 
ter, £l V i S ~ 0.003, and with the amount of bary- 
onic matter, as deduced from primordial nucle- 
osynthesis, fib < 0.05 [0, one concludes that: i) 
most of the matter in the Universe is dark, ii) 
only a small fraction of it is baryonic. 

Particle physics offers a large selection of pos- 
sible candidates for dark matter, once one con- 
siders extensions of the Standard Model (SM). 
The most obvious option is represented by mas- 
sive neutrinos, as suggested by oscillation ef- 
fects observed in solar and atmospheric neutri- 
nos (for the latest data, see Ref.g and Ref.§, 
respectively). In these experiments only differ- 



ences in neutrino squared masses are measured; 
from these results and from the standard formula 
flyh 2 — J2 U m v /(93 eV) a lower bound may be 
established for neutrino contribution to the cos- 
mological matter density, fl v h 2 > 6 • 10 -4 . Obvi- 
ously, the cosmological neutrino contribution may 
be larger, in case of neutrino mass degeneracy. 
However, taking into account the stringent upper 
bound, m Ve < 2.2 eV, derived from Tritium beta 
decay 0, an upper limit on the neutrino cosmo- 
logical contribution may be derived: fl^h 2 < 0.07 
(or fl v < 0.14, using the central value h — 0.7). 

The existence of an upper bound of roughly this 
order of magnitude for the neutrino contribution 
to the matter density in the Universe is not unex- 
pected, since the evolution theory of the primor- 
dial density fluctuations into the present cosmo- 
logical structures indicates that most of the dark 
matter must be comprised of cold particles, i. e. 
of particles that decoupled from the primordial 
plasma when nonrelativistic. Massive neutrinos, 
with masses up to about 1 MeV, are relativis- 
ts at decoupling, and then, by their free stream- 
ing, tend to erase small-scale structures. A re- 
cent investigation || , based on the largest galaxy 
redshift survey (2dFGRS), has derived the upper 
bound fl u < 0.05. 

Thus, by using particle physics results and cos- 
mological observations, one arrives at the conclu- 
sion that light neutrinos are expected to only con- 
stitute a subdominant dark matter component. 

We turn now to cold relics. The order of 
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magnitude of the relic abundance of these cold 
particles, or WIMPs (Weakly Interacting Mas- 
sive Particle), is set by the approximate formula 
Vt x h 2 ~ 3 ■ 10~ 27 cm 3 s -1 / < a a v > int , where 
< a a v >i n t denotes the thermal average of the 
product of the annihilation cross-section a a times 
the relative velocity v for a pair of WIMPs, inte- 
grated from the freeze-out temperature Tf to the 
present-day one |j| . 

As mentioned above, to identify a WIMP with 
a specific particle candidate, one has to move be- 
yond the Standard Model. The only guidance in 
the selection of a specific model is provided by 
theoretical motivations. Indeed, from measure- 
ments at accelerators and in other precision ex- 
periments no significant evidence of physics be- 
yond the SM emerged so far; thus we can only 
count on bounds, which significantly constrain 
the amount of acceptable deviations from the SM 
predictions. 

One of the most strongly motivated extension 
of the Standard Model is represented by Super- 
symmetry. This is the theoretical framework 
within which we consider here possible candidates 
for dark matter. 

2. Supersymmetric models 

Supersymmetric relics are interesting because 
of the following properties: i) if R-parity is 
conserved, the Lightest Supersymmetric Particle 
(LSP) is stable; ii) if colourless and uncharged, 
the LSP is a nice realization of a relic WIMP JToJ . 
In the present paper we assume that supersym- 
mctry exists in Nature and that properties (i) and 
(ii) hold. 

The nature of the LSP depends on the susy- 
breaking mechanism and on the specific regions 
of the susy parameter space considered. Here we 
limit ourselves to gravity-mediated schemes, and 
domains of the parameter space where the LSP is 
the neutralino. The neutralino is defined as the 
lowest-mass linear superposition of photino (7), 
zino (Z) and the two higgsino states (H°, iJf): 
X = 017 + C12Z + CL3H1 + a4-ff|. Hereafter, the 
nature of the neutralino is classified in terms of 
a parameter P, defined as P = a\ + a|. The 
neutralino is called a gaugino when P > 0.9, a 



higgsino when P < 0.1, mixed otherwise. 

Extensive calculations on relic neutralino phe- 
nomenology in gravity-mediated models have 
been performed (among the most recent refer- 
ences on neutralino dark matter see, for instance, 

Even limiting oneself to a the Minimal Su- 
persymmetric extension of the Standard Model 
(MSSM) , a significant variety of different schemes 
exists, from those based on universal or non- 
universal supergravity, with susy parameters de- 
fined at the grand unification scale (GUT), to 
an effective supersymmetric model defined at the 
Electro- Weak (EW) scale. 

2.1. Universal and non universal SUGRA 

The MSSM is described by a Yang-Mills La- 
grangian, by the superpotential, which contains 
all the Yukawa interactions between the stan- 
dard and supersymmetric fields, and by the soft- 
breaking Lagrangian, which models the breaking 
of supersymmetry. The Yukawa interactions are 
described by the parameters h % , which are related 
to the masses of the standard fermions by the 
usual expressions, e.g., nit — h t V2,mi, — h b vi, 
where v\ and v-i are the vev's of the two Higgs 
fields, Hi and H 2 . 

Implementation of this model within a su- 
pergravity scheme leads naturally to a set of 
unification assumptions at a Grand Unification 
(GUT) scale, Mgut- i) Unification of the gaug- 
ino masses: M^Mgut) = ffii/2, ii) Universality 
of the scalar masses with a common mass de- 
noted by mo: ui^Mgut) = mo, iii) Universal- 
ity of the trilinear scalar couplings: A 1 (Mgut) = 
A d (Mgut) = A u (M GU t) = A m . The rele- 
vant parameters of the model at the electro-weak 
(EW) scale are obtained from their corresponding 
values at the Mgut scale by running these down 
according to the renormalization group equations 
(RGE). By requiring that the electroweak sym- 
metry breaking is induced radiatively by the soft 
supersymmetry breaking, one finally reduces the 
model parameters to five: TO1/2, too, Aq, tan/3(= 
V2/V1) and sign u,. In this very strict scheme (uni- 
versal SUGRA) the phenomenology of relic neu- 
tralinos is very sensitive to the way in which vari- 
ous constraints (for instance those on the bottom 
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quark mass, mb, on the top quark mass, m t , and 
on the strong coupling a s ) are implemented. 

Models with unification conditions at the GUT 
scale represent an appealing scenario; however, 
some of the assumptions listed above, particularly 
ii) and iii), are not very solid, because, as was al- 
ready emphasized some time ago [^0| , universality 
might occur at a scale higher than Mqut ~ 10 16 
GeV, e.g., at the Planck scale. 

To take into account the uncertainty in the 
unification scale one may introduce deviations 
in the unification conditions at Mqut- For in- 
stance, deviations from universality in the scalar 
masses at Mqut, which split Mh 1 from Mh 2 
may be parametrized as Mjj.(Mgut) = m o(^ + 
Si). This is the case of non-universal SUGRA 
(nuSUGRA) that we considered in Refs. jL^Q. 
Further extensions of deviations from universality 
in SUGRA models which include squark and/or 
gaiigin o masses are discussed, for instance, in 

The possibility that the initial scale for the 
RGE running, Mj, might be smaller than 
Mqut ~ 10 16 GeV has been raised, on the ba- 
sis of a number of string models (see, for in- 
stance, IjLqpipil and references quoted therein). 
As is stressed in Ref.^U, Mj might be any- 
where between the EW scale and the Planck scale, 
with significant consequences for neutralino phe- 
nomenology, e. g. for the size of the neutralino- 
nucleon cross section. 

2.2. Effective MSSM 

In view of the large uncertainties involved in 
the choice of the scale Mr, the SUGRA schemes 
turn out to be somewhat problematic: the orig- 
inally appealing feature of a universal SUGRA 
with few parameters fails, because of the need to 
take into consideration the variability of Mj or, 
alternatively, to add new parameters which quan- 
tify the various effects of deviation from univer- 
sality at the GUT scale. Thus, it appears more 
convenient to work with a phenomenological susy 
model whose parameters are defined directly at 
the electroweak scale. Here, this effective scheme 
of MSSM is denoted as effMSSM. It provides, at 
the EW scale, a model, defined in terms of a mini- 
mum number of parameters: only those necessary 



to shape the essentials of the theoretical structure 
of MSSM, and of its particle content. Once all ex- 
perimental and theoretical constraints are imple- 
mented in this effMSSM model, one may inves- 
tigate its compatibility with specific theoretical 
schemes at the desired Mj. 

In the effMSSM scheme we consider here, we 
impose a set of assumptions at the electroweak 
scale: a) all trilinear parameters are set to zero 
except those of the third family, which are unified 
to a common value A; b) all squark soft -mass pa- 
rameters are taken degenerate: = m^; c) all 
slepton soft-mass parameters are taken degener- 
ate: mj, = mf, d) the U(l) and 5/7(2) gaugino 
masses, Mi and M2, are assumed to be linked 
by the usual relation Mi = (5/3) tan 2 8\yM2 
(this is the only GUT-induced relation we are 
using, since gaugino mass unification appears 
to be better motivated than scalar masses uni- 
versality). As a consequence, the supersym- 
metric parameter space consists of seven inde- 
pendent parameters. We choose them to be: 
Ma, n, tan 0, itia, rriq, mj, A (tua is the mass of the 
CP-odd neutral Higgs boson). 

Much larger extensions of the supersymmetric 
models could be envisaged: for instance, non- 
unification of the gaugino masses [|l7| , p3| , and 
schemes with CP-violating phases |24|j . 

In the present note we only report results ob- 
tained in the framkework of the effMSSM. For 
a comparison with results derived in other susy 
schemes see Ref.Q. 

2.3. Constraints on supersymmetric pa- 
rameters 

In our exploration of the susy parameter space, 
we have implemented the following experimental 
constraints: accelerators data on supersymmetric 
and Higgs boson searches (CERN e + e~ collider 
LEP2 (2§ and Collider Detector CDF at Fermilab 
pql ); measurements of the b — > s + 7 decay J27]]. 

A further constraint is supplied by the bounds 
on the supersymmetric contribution to muon 
anomalous magnetic moment a™ sy . The Muon 
(g-2) Collaboration at Brookhaven National Lab- 
oratory has recently reported a new experimen- 
tal determination of = (g — 2)/2 p8[ ; this 
confirms their previous measurement ]29| with a 
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precision increased by a factor of two. Unfortu- 
nately, the theoretical evaluation of the SM con- 
tribution to a u , a^ M , is affected by an error which 
is about twice the present experimental uncer- 
tainty. If we use for the leading hadronic contri- 
bution the average of the evaluations in Refs.pof 
|3^| (see also Ref.[j33"f), for higher order correc- 
tions the results of Rcfs.]3^| and for the hadronic 
light-by-light contribution the average of the re- 
sults of Rc fs.|35pq |, we find that a^ M deviates 
from af^ p [[28[ only by a mere 1. 4a, i.e. that 
Aa u ee a e * p - al M sits in the rounded-off 2a 
range: -100 < Aa n • 10 11 < 550. In the present 
paper we restrain a s ™ sy to stay inside this range. 
Parenthetically, we note that Aa u would rise to a 
2.6a level, if for the leading hadronic contribution 
one arbitrarily uses only the results of Ref. p0[ 
and disregards those of Refs.[[3l| [331 ■ 

3. Detection rates for relic neutralinos 

Relic neutralinos may produce many different 
effects, generated either by their interactions with 
ordinary matter or by their pair-annihilations. 
These effects may be detected in experiments de- 
signed to search for WIMPs, cither directly or 
indirectly. 

3.1. Direct detection 

Direct detection aims at the measurement of 
the energy released by a WIMP to a target nu- 
cleus as a consequence of their mutual interac- 
tion p^] . The detection rate is proportional 
to the WIMP-nucleus cross-section convoluted 
over the phase-space distribution function of the 
WIMPs in the galactic halo, calculated at the 
solar-system neighbourhood. As is customary, 
we assume that this distribution function may 
be written as the product of a matter density p x 
times a velocity distribution funcion f(v). Once 
a specific f(v) is chosen, measurements of the 
detection rate in direct WIMP experiments pro- 
vide information on the product of p x times the 
WIMP-nucleus cross-section, or on the product 

(nucleon) i i , . , , . n 

Px "scalar ' wnen coherent interactions domi- 
nate over spin-dependent ones {a scalar denotes 
the WIMP-nucleon scalar cross-section) J3S| . We 
further express p x as p x = £ • pi, where pi is the 



local total dark matter density and £ (£ < 1) is 
a scaling parameter which accounts for the ac- 
tual fraction of local dark matter to be ascribed 
to the candidate \- Thus, in the case of WIMPs 
with coherent interactions with matter, the range 
of WIMP parameters under current exploration 
in WIMP direct searches may approximately be 
given as: 

10~ 10 nbarn < < 2 • lO" 8 nbarn, (1) 

for WIMP masses m x in the range: 

30 GcV <m x < 270 GeV. (2) 

In deriving the ranges in Eqs. a variety 

of WIMP galactic distribution functions and un- 
certainties in astrophysical quantities have been 
taken into account [ |39pC| ]. When interpreted as 
due to a WIMP with coherent interactions, the 
DAMA annual-modulation effect |^l[ provides a 

o ±1, 1 f (nucleon) , . , 

da region m the plane m x — £,a sc , which 
is embedded in the range of Eqs.([j]-||). Other 
experiments of WIMP direct detection provide 

upper bounds on ta^iar^ ©1111] • :t is 
worth stressing that a comparison of results de- 
rived from experiments using different target nu- 
clei is only feasible under the assumption that the 
WIMP-nucleus interaction is dominantly coher- 
ent; otherwise, the comparison is either impossi- 
ble or strongly model-dependent. Interpretation 
of the DAMA annual-modulation results in terms 
of a WIMP with a mixed (coherent and spin de- 
pendent) coupling is discussed in Ref. |Q . 

3.2. Indirect detection 

Indirect signals may be produced by: 

a) Annihilation of neutralinos inside celestial 
bodies (Earth and Sun), where x's have been 
captured and accumulated. Among the vari- 
ous annihilation products, neutrinos would be 
able to emerge form the celestial bodies and be 
detectable as up-going muons in neutrino tele- 
scopes. 

b) Neutralino-neutralino annihilations taking 
place in the galactic halo. Specific signals would 
emerge: most notably, neutrino fluxes, pho- 
ton fluxes, and rare components in cosmic rays 
(positrons, antiprotons, antideuterons). 
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In case of signals of category (a) the detection 
rate is, as in the case of direct detection, pro- 
portional to the product p x 0"^^ on ' (assuming 
that coherent interactions dominate over spin- 
dependent ones); for signals of category (b), the 
detection rates are proportional to the product 
Px < &annV >, where p x is the average of pt{r) 
over the galactic halo. This makes the size of sig- 
nals of type (b) very sensitive to effects of dumpi- 
ness in spatial distribution of WIMPs Jl5| . Actu- 
ally, recent high-resolution simulations of cosmic 
structures appear to favour formation of regions 
of matter overdensity inside galaxies pj| , with 
important implications for neutralino signals (see, 
for instance, |47|] ). 

One may wonder whether constraints on neu- 
tralinos could be derived from some other indi- 
rect means, not related to specific experiments 
designed to search for WIMPs. For instance, it 
has been argued that significant constraints could 
be deduced from helioscismology or from data on 
solar neutrino fluxes Jl8|. However, it has been 
subsequently proved that no constraints can be 
derived at present from solar physics for WIMPs 
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with masses above 30 GeV 49 



4. Results 

A first important question is whether present 
sensitivities in WIMP direct and indirect exper- 
iments are enough to explore some parts of the 
supersymmetric parameter space. Furthermore, 
one wonders whether the relevant neutralinos pro- 
vide a significant contribution to f2 m . These 
two questions are connected, since VL x h 2 and 



nucLeon 



scalar 



(and consequently also the neutralino 
direct detection rate) are approximately anticor- 
related |o||5l). 

Fig.l provides an answer to the above ques- 
tions. First, it shows that present experiments 
of WIMP direct search actually explore signif- 
icant regions of the supersymmetric parameter 
space. This property was already pointed out in 
Ref. ^2j; actually, in the mid-nineties, with the 
advent of new detectors, WIMP direct searches 
acquired the capability of exploring regions of 
the susy parameter space not yet excluded by ac- 
celerator bounds. Since then, the progressively 




Figure 1. Scatter plot of cr^"" ; ^°™ J versus fl x 
in the effective MSSM model described in the 
text. Only configurations with positive p and 
with m x in the range of Eq.(Q) are shown. The 
two horizontal lines bracket the sensitivity region 
for the WIMP-nucleon scalar cross-section de- 
fined in Eq. (|l|) . The two vertical lines denote the 
range 0.05 < Sl x < 0.3. Dots denote gauginos, 
circles denote higgsinos and crosses denote mixed 
configurations. 



more stringent limits established by new acceler- 
ator data have been counterbalanced by the im- 
proved sensitivities of WIMP direct experiments 
|Tl-I3|. For analyses by other groups see for in- 
stance Refs.(l|-[l9|. 

Secondly, Fig.l shows that exploration by 
WIMP direct search is effective also in sectors 
where relic neutralinos may provide a substantial 
contribution to Q m . In view of the anticorrela- 
tion between fl x h 2 and cr["^^ on \ this property 
is far from being trivial. 

However, also the case of relic neutralinos 
which only constitute a subdominant dark mat- 
ter population are interesting in view of their ex- 
plorability by WIMP direct experiments ||50| , pl| . 
To further illustrate this point, in Fig. 2 we dis- 
play a scatter plot of the susy configurations 
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Figure 2. Scatter plot of p x versus D, x h 2 . This 
plot is derived from the value [p x /(0.3 GeV cm -3 ) 

' CT sc^ar = 1 • 10~ 9 nbarn and by taking m x 

in the range of Eq. (||) , according to the procedure 
outlined in the text, in case of effMSSM. The two 
horizontal lines delimit the range 0.2 GeV cm~ 3 
< Px - I-? GeV cm~ 3 Q. The two vertical lines 
delimit the range 0.05 < Q x < 0.3. Dots denote 
gauginos, circles denote higgsinos and crosses de- 
note mixed configurations. 



which belong to the value [p x /(0.3 GeV cm 3 ) 



c^ar° ]expt =1-10 9 nbarn, which is repre- 
sentative of a value of p x ■ cr^™^ 01 ^ inside the 
experimental sensitivity range of Eq. ([!]), and to 
the mass interval of Eq.(|J). The plot is ob- 
tained in the following way: 1) p x is evaluated 

r (nuclcon) I / (nucleoli) t (nuclcon) . t 

aS [Px ' ^scalar \expt / Scalar ^scalar IS Cal- 
culated within the susy model), 2) to each value 
of p x one associates the corresponding calculated 
value of il x h 2 . For a full discussion of the var- 
ious sectors displayed in Fig. 2 we refer to Jsi] . 
Here we wish to call the attention only to two 
of these sectors. The most interesting one is 
obviously the rectangular (cross-hatched) region, 
where 0.2 GeV cm" 3 < p x < 1.7 GeV cm" 3 Q 
and 0.05 < f2 m /i 2 < 0.3; this corresponds to a sit- 



uation where the neutralino may be a dominant 
dark matter component. A second interesting re- 
gion is the band delimited by the slanted dot- 
dashed lines and simply-hatched in the figure. 
For configurations which fall inside this band, the 
neutralino would provide only a fraction of the 
cold dark matter, at the level of local density and 
of the average relic abundance; a situation which 
would be possible, for instance, if the neutralino 
is not the unique cold dark matter particle com- 
ponent. To neutralinos belonging to these config- 
urations one should assign a rescaled local density. 
Our results prove that also these relic neutralinos 
are actually probed by WIMP direct searches. 

We turn now to indirect WIMP searches. It 
may be shown that also measurements of up- 
going muons at neutrino telescopes ]53|-|5^] might 
be capable of exploring interesting neutralino 
configurations ]57[ | , with some overlap with those 
analysed in WIMP direct detection experiments 
f39| . Fig. 3 displays the scatter plot for the up- 
going muon flux from the center of the Earth, 
evaluated in effMSSM, for a standard galactic dis- 
tribution function and taking into account the ef- 
fect of — > v T oscillations, for neutrino mass 
and mixing angle parameters which best fit the 
atmospheric neutrino data fH5cf|. 

A number of susy configurations are at the level 
of current experimental bounds. However, one 
should be aware of the fact that the efficiency of 
capture of WIMPs by the Earth, and then also 
the up-going muon flux from the center of the 
Earth, are very sensitive to the modifications in- 
duced by the Sun on the local low-energy WIMP 
population. Depending on the Damour-Krauss 
conjecture |^| or the Gould- Alam conjecture Q), 
the up-going muon flux from the center of the 
Earth might be either enhanced or suppressed, 
as compared to the standard one. Therefore, at 
present, derivations of constraints on supersym- 
mctric configurations from experimental bounds 
on up-going muon flux from the center of the 
Earth is very hazardous. This important point 
is often disregarded in the literature. 

Very promising are also strategies for mea- 
suring in space signals due to neutralino pair- 
annihilations taking place in the halo. Here we 
briefly mention only some signals consisting of 
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Figure 3. Scatter plot for the up-going muon flux 
due to neutralino-neutralino annihilation in the 
center of the Earth in effMSSM. Dots (crosses) 
denote configurations with Q x h 2 < 0.05 (0.05 < 
fl x h 2 < 0.3). The dashed line denotes the 90 
% C.L. upper bound of Rcf . |54j] , the solid line 
denotes the 90 % C.L. upper bound of Ref . |55|| . 



rare components in cosmic rays. 

One of the most widely discussed case is a pos- 
sible primary flux of low-energy antiprotons. Cal- 
culations of secondary antiprotons, due to inter- 
actions of cosmic rays with interstellar material, 
have greatly improved in recent years J62| , |63[ and 
have shown that at low-energy (antiproton kinetic 
energy T p less then about 1 GeV) the secondary 
spectrum is much flatter than previously believed 
and fits remarkably well the experimental data 
64 1 . Therefore, only tiny room is left to exotic 
components, and the discrimination between sec- 
ondary and primary components becomes very 
problematic. However, since a number of super- 
symmetric configurations reach or even exceed 
the level of the experimental spectrum 62 1, high 
statistics in p measurements and great accuracy 
in the evaluation of the primary and secondary 
spectra may provide useful constraints for super- 
symmetry. 

Recently, it has been suggested that a measure- 
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Figure 4. Top-of-atmosphere antideuteron en- 
ergy spectra from Ref. pq| . The solid line dis- 
plays the secondary flux, the other lines denote 
primary fluxes for four different neutralino com- 
positions (see Ref. |6j| for details). 



mcnt of low-energy antideuterons in space would 
provide a way of looking at a distinctive signa- 
ture of WIMP pair annihilation in the halo Q . 
In fact at energies below about 1 GeV per nucleon 
the primary spectrum would be quite dominant 
over the secondary one (see Fig. 4). This prop- 
erty is at variance with what happens in case of 
the antiproton low-energy spectra. Thus, the an- 
tideuteron measurements open up quite new in- 
teresting perspectives for detecting x — X annihi- 
lation in space. 

We wish also to mention that the Heat Col- 
laboration |^] reported evidence for an excess 
of positrons in cosmic rays at energies around 
10 GeV. This might be considered as due to 
neutralino-neutralino annihilation (see for in- 
stance |67]]). However, it should be noted that 
this interpretation requires an ad hoc boost fac- 
tor in order to fit the experimental data. 

Other important signals possibly produced by 
neutralino pair-annihilations in the halo are pro- 
vided by gamma-rays; for this class of processes 
see, for instance, Rcfs.[b8pil. 
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In conclusion, it may be said that quite differ- 
ent strategies of searching for WIMPs are now 
capable of providing independent information on 
supersymmetric relic particles. 
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